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The invention concerns a method for focusing acoustic waves useful for obtaining an image of a field to 
be observed in a dissipative heterogeneous medium (2, 3) around which acoustic transducers (T1-Tn, 
n-T'm) forming an imaging network and a target network. The method consists in following a training 
step during which pulse responses from the medium are measured between each transducer (Ti) of the 
imaging network (5) and several transducers (Tj) of the target network (6); deducing therefrom reference 
signals to be emitted by the transducers of the imaging network to produce a focused acoustic pulse in 
each transducer of the target network, then cumulatively, in determining reference signals to be emitted 
to focus an acoustic pulse on predetermined points in the medium. Said reference signals are stored 
and used subsequently to generate an acoustic image of the medium. 
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PROCEDE ET DISPOSITIF NON INVASIF DE FOCALISATION D'ONDES ACOUSTIQUES. 

Proc^d^ de focalisation d'ondes acoustiques utilisable 
pour obtenir une Image d*un champ k observer dans un mi- 
lieu h6t6rogene disslpatif (2, 3) autour duquei on fixe des 
transducteurs acoustiques (Tl-Tn, Tl-Tm) formant un r6- 
seau d'imagerie et un r§seau clble. On suit une 6tape d'ap- 
prentlssage au cours de laquelle on mesure des reponses 
Impulsionnelles du milieu entre chaque transducteur (Ti) du 
reseau d'imagerie (5) et plusieurs transducteurs (Tj) du r6- 
seau cible (6). On en deduit des signaux de r6f6rence ^ 
emettre par les transducteurs du r6seau d'imagerie pour 
produire une impulsion acoustique focalis6e en chaque 
transducteur du r6seau cible, puis de proche en proche, on 
determine des signaux de r6f6rence a 6mettre pour focali- 
ser une impulsion acoustique sur des points predetermines 
dans le milieu. Ces signaux de reference sont memorises et 
utilise ensuite pour g^n^rer une Image acoustique du mi- 
lieu. 
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Procede et dispositif non invasif de f ocalisa1:ion ondes 
acous-tigues . 

La presente invention est relative aux precedes et 
5 dispositif s non invasif s de focalisation d' ondes 
acoustiques, notanmient ultrasonores . 

Plus particulierement , 1' invention concerne un 
procede non invasif de focalisation d' ondes acoustiques 
dans un milieu het6rog6ne dissipatif comprenant un milieu 

10 sensiblement homogSne (par exemple, le cerveau) entour6 au 
moins partiellement par une couche aberratrice dissipative 
(par exemple, le crSne) qui genere des aberrations dans la 
propagation des ondes acoustiques, les ondes acoustiques 
etant 6mises depuis I'exterieur de la couche aberratrice et 

15 focalis6es dans le milieu sensiblement homogfene. 

Les precedes de ce type qui sont courammeht 
utilises ne permettent pas d'obtenir une bonne focalisation 
des ondes acoustiques a I'interieur du milieu, et lorsque 
ces precedes sont utilises dans des applications 

20 d'imagerie/ ils ne permettent done pas d'obtenir une bonne 
resolution et un bon contraste d' image lorsque les 
aberrations de propagation sont import antes, par exemple 
lorsqu'on realise une echographie du cerveau depuis 
l'ext§rieur du crane. 

25 La presente invention a notamment pour but de 

pallier cet inconvenient. 

A cet effet, selon 1' invention, un procede de 
focalisation du genre en question est caract6ris6 en ce 
qu'il comporte les etapes suivantes : 

30 (a) une etape initiale de positionnement au cours 

de laquelle on fixe un nombre t superieur ^ 2 de 
transducteurs acoustiques dans des positions pr6d6termin6es 
a I'exterieur de la couche aberratrice, ces transducteurs 
etant en contact (direct ou indirect) avec ladite couche 

35 aberratrice et formant au moins : 



2815717 



2 

- un reseau d' imager ie qui regroupe un nombre n 
compris entre 1 et t desdits transducteurs, 

- et un reseau cible qui regroupe un nombre m 
compris entre 1 et t desdits transducteurs (ces deux 
reseaux peuvent etre enti^rement distincts, ou comporter 
certains transducteurs communs, ou encore comporter chacun 
la totalite des transducteurs susmentionnes) , 

(b) une 6tape d'apprentissage comprenant elle-m§me 
les sous 6tapes suivantes : 

(bl) une sous-6tape d' apprentissage de la 
focalisation du reseau d' imagerie sur le reseau cible, 
sous-6tape au cours de laquelle : 

(bll) on determine des r6ponses 

impulsionnelles hri (t) du milieu h6t6rog6ne dissipatif, 
respectivement entre chaque transducteur i du r6seau 
d' imagerie et plusieurs points de focalisation r situes sur 
la couche aberratrice en correspondance respective avec des 
transducteurs du reseau cible (cette determination peut 
etre faite par mesure directe si I'on fait 6raettre des 
impulsions acoustiques par les transducteurs du reseau 
cible, ou le cas 6ch6ant par mesure et calcul si I'on fait 
emettre des signaux acoustiques autres que des impulsions 
par les transducteurs du reseau cible, les valeurs ainsi 
mesurees et/ou calcul^es pouvant ensuite etre le cas 
ech^ant corrigees par repropagation num^rique pour simuler 
des transducteurs situ6s directement au contact de la 
couche aberratrice si les transducteurs ne sont pas au 
contact direct de la couche aberratrice), ces r6ponses 
impulsionnelles etant m6moris6es sous forme numerique avec 
un certain echantillonnage temporel qui determine un nombre 
P de composantes f requentielles de la r^ponse 
impulsionnelle, de frequences respectives tok, i 6tant un 
indice compris entre 1 et n qui designe un transducteur du 
reseau d' imagerie, r etant un indice compris entre 1 et m 
qui designe un point de focalisation corrfispondant A un 
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transducteur du reseau cible et k etant un indice compris 
entre 1 et p qui design© une composante f requentielle, 

(bl2) a partir de ces reponses 
impulsionnelles, on calcule, pour chaque point de 
focalisation r correspondant a un transducteur du r6seau 
cible ^ un ensemble de n signaux de reference temporels 
e'i(t,r), i variant entre 1 et n, tels que, si la parol 
aberratrice etait enlevee au voisinage du point de 
focalisation r, 1' emission de ces signaux de reference par 
les diff^rents transducteurs i du reseau d' imagerie 
g6n6rerait un signal predetermine (par exemple une 
impulsion acoustique) focalise sur le point de focalisation 
r, 

(b2) une sous-etape de focalisation en un nombre 
R de points de focalisation pr6d6termin6s situ6s dans le 
milieu sensiblement homogene, d* indices q compris entre m+1 
et m+R, cette sous-6tape consistant a determiner pour 
chacun de ces points de focalisation q, en s'eloignant pas 
a pas des points de focalisation lam correspondant aux 
transducteurs du reseau cible, des signaux de r6f6rence 
e'i(t,q) a faire 6mettre par les differents transducteurs i 
du reseau d'imagerie pour gen6rer une impulsion focalisee 
sur ledit point de focalisation q, les signaux de r6f6rence 
e'i(t,q) etant determines pour chaque point de focalisation 
q en procedant comme suit : 

(b21) une premiere estimation de e'i(t,q), 
pour i allant de 1 a n, est calculee a partir d'au moins un 
signal de reference e'i(t,qO), qO etant 1 ' indice d'au moins 
un point de focalisation proche du point de focalisation q 
pour lequel le signal de reference a deja ete determine, ce 
calcul etant fait en utilisant une celerite moyenne des 
ondes acoustiques dans le milieu sensiblement homogene, 

(b22) on fait emettre par les transducteurs 
du reseau d'imagerie, par iterations, les estimations 
precedemment obtenues des signaux de reference e'i(t,q). 
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puis on capte avec les m^mes transducteurs des signaux 
Si(t,q) r6trodiffus6s par le milieu het6rogene dissipatif, 
puis on modifie pour 1' iteration suivante ces signaux de 
reference e'i(t,q) de la maniere suivante : 

ou les valeurs a,(q) et x,iq) sont un facteur d' amplitude 
et un retard correctifs, calcules pour maximiser un crit^re 
de coherence c entre lesdits signaux r6trodif fus6s, 
lesdites iterations 6tant arrgt^es lorsque le critdre C 
attaint un seuil predetermine, 

^^^^ °" memorise les signaux de reference 
e'i(t,q), au moins pour q compris entre m+1 et m+R, 

(c) et une etape de focalisation au cours de 
laquelle, pour au moins un desdits points de focalisation 
q, on fait emettre respectivement par les transducteurs du 
reseau d'imagerie, lesdits signaux de reference e'i(t,q), i 
6tant un indice compris entre 1 et n designant un 
transducteur du reseau d'imagerie. 

Grace a ces dispositions, on s'affranchit des 
aberrations de propagation des ondes acoustiques dans le 
milieu heterogene dissipatif, et on obtient une 
focalisation tr^s precise qui peut notamment permettre 
d'obtenir une image echographique fidele et precise d'un 
champ a observer A travers la couche aberratrice par 
25 retrodif fusion, lorsqu'on emet successivement des ondes 
acoustiques focalisees sur differents points du champ ^ 
observer et que I'on capte les ondes acoustiques 
retrodif fusees . 

Cette focalisation precise peut egalement etre 
utilisee dans d'autres applications que 1' echographie, 
notamment : 

I'imagerie doppler couleur, 

les methodes d'imagerie par eiastographie, 
telles que celle d6crite dans le document WO-A-00/55 616, 

les methodes d'imagerie non lineaire ("harmonic 



20 
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35 



INSDOaO: <Fn ai5717A1_l_> 



2815717 



5 

imaging") , 

les m^thodes de traitement par destruction 
localisee d'une partie du milieu heterogene dissipatif, 
notamment par hyperthermie, 
5 - les methodes de mesure de paramfetres 

d' absorption optique des tissus avec activation par 
ultrasons, etc. 

Dans des modes de realisation pref6r6s de 
1' invention, on peut 6ventuellement avoir recours en outre 

10 a I'une et/ou h 1' autre des dispositions suivantes : 

au cours de la sous-6tape (bll), lorsqu' au moins 
certains transducteurs (du r6seau cible et/ou du r6seau 
d'imagerie) sont en contact avec un milieu homogdne 
interm6diaire (par exemple un gel) lui-mSme en contact avec 

15 la couche aberratrice, on corrige les reponses 
impulsionnelles hri(t) par repropagation numerique pour 
simuler des transducteurs situes directement au contact de 
la couche aberratrice ; 

la sous-etape (bl2) comporte elle-meme les sous- 

20 etapes suivantes : 

(bl21) on determine p matrices de transfert 
H(cok) = [Hri (CDk) ] , i allant de 1 a n et r allant de 1 a m, 
ou Hri (ok) est la valeur, a la frequence ©k, de la 
transformee de Fourier de la reponse impulsionnelle hri(t), 

25 (bl22) on determine pour chaque point de 

focalisation r correspondant a un transducteur du reseau 
cible, n composantes Ei(G)k,r), i variant entre 1 et n, 
telles que F(©k, r) =H (CDk) .E (cok, r) , ou E((ok,r) = [Ei(cok,r)] 
est un vecteur A n composantes, F(a)k,r) est un vecteur a m 

30 composantes Fl((Dk,r), 1 variant entre 1 et m, ces m 
composantes Fl(a)k,r) correspondant a une focalisation 
souhaitee des ondes acoustiques a la frequence cok sur le 
point de focalisation r correspondant a un transducteur du 
reseau cible. 
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(bl23) on en d^duit, pour chaque point de 
focalisation r correspondant a un transducteur du r6seau 
cible, un vecteur de n signaux temporels e (t, r) = [ei (t , r ) ] , 

i variant entre 1 et n, o£x e:M=±Ei(ojk,r).e"^' en notation 

complexe, ces signaux ei(t,r) 6tant adap^es pour que 
1' Emission de ces derniers respectivement par les 
diff6rents transducteurs i du r6seau d' imager ie gen^re une 
impulsion acoustique focalis6e sur le point de focalisation 
r du r6seau cible, 

(bl24) une sous-6tape de correction des 
aberrations g6n6r6es par la couche aberratrice entre le 
milieu sensiblement homogfene et chaque transducteur cible 
r, ces aberrations 6tant estim^es sur la base des mesures 
pr6c6demment effectu6es, les aberrations ainsi estim6es 
6tant utilis6es pour calculer lesdits signaux temporels de 
reference e'i(t,r) ; 

- au cours de la sous-6tape (bl22) on calcule p 
matrices H-^(G)k), respectivement par r6gularisation et 
inversion des matrices de transfert H{<ok) , et pour chaque 
transducteur r du r6seau cible, on calcule le vecteur 
E(Q)k, r) par la formule : 

E(cok,r) =H-Ma)k) .F(a)k, j) ; 

- au cours de I'^tape (bl22) , les" composantes 
Fl«Dk,r) du vecteur F((ok,r) correspondant a la repartition 
spatiale du champ desire a la frequence tok, sont 6gales ^ 
0 pour l?tr et egale a 1 pour l=r ; 

- au cours de la sous-6tape {bl24), on assimile la 
parol aberratrice, au voisinage de chaque point de 
focalisation r correspondant a un transducteur du r6seau 
cible, a un filtre ^ reponse impulsionnelle finie, d6fini A 
chaque frequence <ok par une amplitude Gr(a)k) et une phase 
<^r«ok), la sous-6tape (bl24) comportant elle-meme les sous- 
etapes suivantes : 
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(bl241) on calcule, pour chaque frequence <ok, 
1' amplitude Gr((Dk) et la phase <|)r(<Dk), a partir soit des 
signaux ei(t,r), soit des vecteurs E((Ok, r) , 

(bl242) on calcule p matrices de transfert 
5 corrigees H' (o>k) = [H' j i (ook) ] , ou 

(bl243) on determine pour chaque transducteur r 
du reseau cible^ n composantes E^Kcok^r), i variant entre 
1 et n^ telles que F (<ok, r) =H' (©k) . E' (©k, r) , oil E'(G>k,r) = 

10 [Ei(<Dk,r)] est un vecteur A n composantes, F(cok,r) est un 
vecteur ^ m composantes Fl(cok,r), 1 variant entre 1 et m, 
ces m composantes Fl(a>k,r) correspondant a une 
focalisation souhait6e des ondes acoustiques a la frequence 
cok sur le point de focalisation r correspondant a un 

15 transducteur du reseau cible, 

(bl244) on en deduit, pour chaque point de 
focalisation r correspondant ^ un transducteur du r6seau 
cible, un vecteur de n signaux temporels de reference 
e' (t, r) = [e' i (t^r) ] , i variant entre 1 et n, oCi 

20 e'i(t,r)=^£Pi{a^,r)^^'^' en notation complexe ; 

au cours de la sous-6tape (bl241) , on calcule 
1' amplitude Gr((ak) et la phase <|)r (©k) corame suit : 



Ei(a}k^)Mi (ok^) 




2 5 ^H^O=ig(arg(£(a>*,K)))-ar^^^^ 



ou 



Ei* est la valeur complexe conjuguee de Ei, 

et AT(rO,r,i)=(d(rO,i)-d(r,i))/c, d(r,i) etant la 
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distance entre le transducteur i et le point de 
focalisation r, et d(rO,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO ; 

la sous-etape (bl2) comporte elle-meme les sous- 
etapes suivantes : 

(bl21) on determine p matrices de transfert 
H«Dk) = [Hri(o)k)], i allant de 1 A n et r allant de 1 a m, 
ou Hri((Dk) est la valeur, a la frequence ok, de la 
transformee de Fourier de la reponse impulsionnelle hri(t), 

(bl22') on corrige les matrices de transfert 
H(cok) pour s'affranchir des aberrations g6ner6es par la 
pairoi aberratrice au voisinage de chaque point de 
focalisation r, cette correction 6tant effectuee a partir 
des r6ponses impulsionnelles hri(t) prec6denmient 
determinees, et on obtient ainsi des matrices de transfert 
corrigees H' (oak) , 

(bl23' ) on determine pour chaque point de 
focalisation r correspondant A un transducteur du r6seau 
cible, n composantes E'i(CDk,r), i variant entre 1 et n, 
telles que F(a>k,r)=H^ (cok) .E' (cok, r) , ou E' (<ok,r) = 
[E'i(cok, r)] est un vecteur a n composantes, F((ok,r) est un 
vecteur a m composantes Fl(©k,r), 1 variant entre 1 et m, 
ces m composantes Fl (cok, r) correspondant a une 
focalisation souhaitee des ondes acoustiques a la frequence 
cok sur le point de focalisation r correspondant ^ un 
transducteur du reseau cible, 

(bl24') on en d6duit, pour chaque point de 
focalisation r correspondant a un transducteur du r6seau 
cible, un vecteur de n signaux teinporels 
e' (t,r)=[e'i(t,r) ] , i variant entre 1 et n, ou 

e'i{t,r)=^Ei{a}k,r)^^'^' en notation complexe, les signaux 

e'i(t,r) etant lesdits signaux de reference ; 

- au cours de la sous-etape (bl23' ) on calcule p 
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matrices H'-^(o)k), respect ivement par regularisation et 
inversion des matrices de transfert H' (wk) , et pour chaque 
transducteur r du reseau cible, on calcule le vecteur 
E' {<Dk,r) par la formule : 

E' (ok, r ) =H' (ottk) . F (ok, j ) ; 

- au cours de I'etape (bl23'), les composantes 
Fl(ok,r) du vecteur F(ok,r) correspondant a la repartition 
spatiale du champ desire a la frequence ok, sont 6gales a 
0 pour l?tr et egale A 1 pour l=r ; 

- au cours de la sous-etape (bl22' ) , on assimile 
la parol aberratrice, au voisinage de chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, A un filtre ^ reponse impulsionnelle finie, defini k 
chaque frequence ok par une amplitude Gr(ok) et une phase 
«j>r(cok), la sous-6tape (bl22' ) comportant elle-m§me les 
sous-6tapes suivantes : 

(bl22'l) on calcule, pour chaque frequence 
ok, 1' amplitude Gr (ok) et la phase <|.r(ok), a partir des 
reponses impulsionnelles pr6c6demment d6termin6es, 

(bl22'2) on calcule p matrices de transfert 
corrigdes H' (mk) = [H' ji (ok) ] , oh 

- au cours de la sous-4tape (bl22'l) on calcule, 
pour chaque frequence ok, 1' amplitude Gr(ok) et la phase 
<|>c(ok), de la maniSre suivante : 



COk)Jfrt)J (Wk) 



2815717 



10 

H*ri designe la valeur complexe conjugu6e de Hri, 
et AT(rO,r,i)=(d(rO,i)-d(r,i))/c, d(r,i) 6tant la 
distance entre le transducteur i et le point de 
focalisation r, et d(rO,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO ; 

- au cours de l'6tape (c) , on fait suivre la sous- 
etape (cl) par les sous-§tapes suivantes : 

{c2) on fait capter par lesdits transducteurs du 
reseau d'imagerie des signaux Si(t) retrodif fus6s par le 
milieu h6t6rog6ne dissipatif, 

(c3) on convolue le signal de r6f6rence 6mis par 
cheque transducteur du reseau d'imagerie avec le signal 
retrodiffus6 capt6 par ce transducteur, 

(c4) puis on somme les produits de convolution 
ainsi obtenus, 

1' etape (c) §tant renouvel^e pour une plurality de points 
situ6s dans le milieu sensiblement homogdne ; 

- au cours de la sous-6tape (b21) , la premiere 
estimation de chaque signal de r6f6rence est e'i(t,q)= 
e'i(ts+ei(q) ,qO) pour chaque point de focalisation q, qQ' 
etant I'indice d'un point de focalisation proche du point q 
pour lequel le signal de reference a deja 6te d6termin6, 
ei(q) etant un retard 6gal A une valeur 5i(q)/c, oCi c est 
la c616rite moyenne des ondes acoustiques dans le milieu, 
et 8i(q) est 6gal A une difference entre d'une part, une 
distance entre le transducteur i du reseau d'imagerie et le 
point de focalisation qO, et d' autre part, une distance 
entre le transducteur i du r6seau d'imagerie et le point de 
focalisation q, 

- au cours de la sous-6tape (b2), lorsqu'au moins 
certains transducteurs d' indice v du reseau d'imagerie ne 
sont pas directement au contact de la couche aberratrice, 
on corrige les signaux e'Wt,q) correspondants par 
repropagation numerique pour simuler des transducteurs 
places en contact direct avec la couche aberratrice ; 



I 
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- au cours de la sous-etape (b22), on recherche 
les valeurs ai(q) et Ti(q) pour raaximiser le crit6re de 
coherence C suivant : 



C = 





n 


2 


< 




> 



ou : 



5 • ^i(^^) = •y,(0®^i(^^^) f ® representant 1' operation de 

convolution, 

et <> represente une moyenne temporelle ; 

- au cours de la sous-etape (b22) , les valeurs 
tiCq) sont calculees en maximisant une fonction 

10 d' intercorrelation, pour des transducteurs voisins du 
reseau d' imager ie, des signaux gi(t,q) et gi+i(t,q) ; 

" au cours de la sous-etape (b22) , les valeurs 
oti(q) sont calcul§es de maniere a egaliser sur I'indice i 
1' amplitude maximale des fonctions gi(t,q) ; 

15 - au cours de la sous-etape (b22), les valeurs 

ai(q) et Xi(q) sont calculees en realisant une 
intercorr61ation, pour des transducteurs voisins du reseau 
d'imagerie, des signaux gi(t,q) et gi+i(t,q) ; 

- au cours de la sous-etape (b22), les valeurs 
20 ttiCq) et Ti(q) sont calculees de maniere a egaliser sur 

I'indice i 1' amplitude maximale des fonctions gi(t,q) ; 

la s6us-6tape (b22) relative a chaque point de 
focalisation q est r6alisee immediatement aprfes la sous- 
etape {b21) relative au meme point de focalisation q ; 

25 - le milieu heterogene dissipatif est constitue 

par le cerveau entoure par le crane ; 

le reseau d' imager ie et le reseau cible sont 
deux reseaux distinct s disposes de part et d' autre du 
milieu heterogene dissipatif ; 

30 - tous les transducteurs appartiennent a la fois 

au reseau d' imagerie et au reseau cible ; 
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les ondes acoustiques sont des ondes 
ultrasonores. 

Par ailleurs, 1' invention a 6galement pour objet un 
dispositif adapts pour mettre en oeuvre le precede d^fini 
ci-dessus . 

D'autres caract6ristiques et avantages de 
!• invention apparaitront au cours de la description 
suivante d'une de ses formes de realisation, donn6e k titre 
d'exemple non limitatif, en regard du dessin joint. 

Sur le dessin, la figure unique represente un 
dispositif d'imagerie ultrasonore selon une forme de 
realisation de 1' invention. 

Le dispositif 1 d'imagerie ultrasonore represents 
sur le dessin est adapte pour r6aliser une image 
echographique par ultrasons du cerveau 2 d'un patient (A 
des frequences par exemple de I'ordre de 1 A 3 MHz), depuis 
I'exterieur du crane 3, le cerveau 2 constituant un milieu 
sensiblement homogene pour la propagation des ondes 
acoustiques et le crSne 3 constituant une couche 
aberratrice dissipative, de sorte que 1' ensemble de la 
bolte crSnienne 2,3 constitue un milieu heterogSne 
dissipatif . 

En variante, 1' invention serait applicable 
notamment : 

a 1 • imagerie acoustique de tout autre milieu 
hetirogene dissipatif non homogene comprenant un milieu 
sensiblement homogene entoure par une couche dissipative 
relativement peu epaisse g^nerant des aberrations dans la 
propagation des ondes ultrasonores, 

ou i tout autre precede impliquant au moins une 
focalisation a 1' emission dans un tel milieu. 

Dans 1' exemple represente sur le dessin, le 
dispositif d'imagerie 1 comporte un micro-ordinateur 4, ou 
tout autre dispositif de commande et/ou de visualisation 
des images ultrasonores, ce micro-ordinateur comportant 
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classiquement un clavier 4a associe eventuellement a 
d'autres interfaces de coinmande et un ecran 4b permettant 
de visualiser les images du cerveau 2. 

Par ailleurs, le dispositif d'imagerie 1 comporte 
deux reseaux 5, 6 de transducteurs ultrasons Tl, T2.„ Tn et 
T'l, T'2... T'ra formant par exemple deux barrettes lineaires 
de transducteurs qui sont disposees de part et d' autre du 
crane 3 de 1 ' utilisateur, dans des positions g6ometriques 
predeterminees I'une par rapport a 1' autre, chaque barrette 
de transducteur 5, 6 etant mise en contact avec le crane 3 
par 1 • intermediaire d'une couche 7 de gel ou similaire. 

Les differents transducteurs Tl, T2... Tn et T"l, 
T'2... T'm peuvent etre commandes directement par le micro- 
ordinateur 4, ou de preference par une unit6 centrale 
electronique CPU contenue par exemple dans une baie 
electronique 8 et elle-m§me commandee par le micro- 
ordinateur 4 . 

Avantageusement, chacun des transducteurs Tl, T2... 
Tn, T'l, T'2, T'm est relie a un echantillonneur, 
respectivement El, E2... En, E'l, E'2, E'm, et chaque 
echantillonneur est lui-meme relie a une m6moire, 
respectivement Ml, M2... Mm, M'l, M'2... M'm et a une unit6 
centrale CI, C2, ...Cm, C 1 , C 2 , ... Cm. Ces memoires et ces 
unites centrales sont a leur tour relies, directement ou 
indirectement , a 1' unite centrale CPU susmentionnee^ qui 
est par ailleurs reliee au moins a une memoire centrale M. 

Le dispositif qui vient d'etre decrit fonctionne 
comme suit . 

Initialement , les deux reseaux de transducteurs 5, 
6 sont fixes de part et d' autre du crane 3 du patient, dans 
lesdites positions predeterminees. A cet effet, les reseaux 
de transducteurs 5, 6, dits respectivement reseau 
d'imagerie et reseau cible, peuvent etre portes par un 
support rigide tel qu'un casque (non represente) dispose 
autour de la tete du patient. 



2815717 



14 

Ensuite, le dispositif suit une 6tape 
d'apprentissage de quelques minutes (avantageusement 1 ^ 3 
minutes) permettant de prendre en compte I'ensemble des 
aberrations de propagation dues au caractdre non homog^ne 
du milieu dissipatif forme par le crane 3 et le cerveau 2. 

Au cours de cette 6tape d'apprentissage on fait 
d'abord 4mettre successivement par chacun des transducteurs 
Tl, T2... Ti,... Tn du r^seau d'imagerie 5, une impulsion 
acoustique, et pour chaque impulsion 6mise par 1-un des 
transducteurs Ti du r^seau d'imagerie, on enregistre le 
signal capte par les transducteurs T'l, T'l... T'r,.„ T'm du 
reseau cible 6, c'est-a-dire la r^ponse impulsionnelie 
hri(t) du milieu h6t6rog§ne dissipatif entre le 
transducteur i consid6r6 du reseau d'imagerie 5 et chaque 
transducteur j du reseau cible 6. 

Chaque r^ponse impulsionnelie hri(t) est 
enregistree sous forme num^rique avec un certain- 
echantillonnage temporel qui determine un certain nombre p 
de composantes f requentielles monochromatiques de la 
r^ponse impulsionnelie, correspondant chacune A une 
frequence <ok, k etant un indice compris entre 1 et p. 

Dans le cas envisage ici, oia au moins certains 
transducteurs du reseau cible et/ou du reseau d'imagerie ne 
sont pas directement au contact avec la couche aberratrice 
3, on corrige les reponses impulsionnelles pour simuler des 
transducteurs virtuels disposes au contact de ladite couche 
aberratrice. La position de la couche par rapport aux 
transducteurs peut 6ventuellement etre obtenue par imagerie 
conventionnelle (echographie ultrasonore, scanner X, IRM, 
etc.). Les r6ponses impulsionnelles corrigees sont 
calcul6es par un algorithms de repropagation num^rique 
connu, d^crit notamment dans les articles suivants : 

"Ultrasonic beam steering through inhomogeneous 
layers with a time reversal mirror", C.DORME, M. FINK, IEEE 
Transactions Ultrasonics^ Ferroelectric and Frequency 
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Control, 43 (1), janvier 1996, p 167-175, 

"Focusing and steering through absorbing and 
aberrating layers : Application to ultrasonic propagation 
through the skull" Journal of Acoustical , Society of 
America, 103 (5), Mai 1998, p. 2403-2410, 

et "Propagation and backpropagation for 
ultrasonic wavefront design" Liu, D.-L., and Waag, R- C. 
IEEE Trans, on Ultras. Ferro. and Freq. Contr. 44(1): 1-13 
(1997). 

Dans ce qui suit, hri(t) denoiranera done les reponses 
impulsionnelles pour des 616ments (reels ou virtuels) 
situ6s centre la couche aberratrice. De plus, les 616ments 
virtuels ou r6els situ6s centre la couche aberratrice 3 
seront appel6s ci-apr6s "points de f ocalisation" d' indice r 
compris entre 1 et m, 

Lorsque les transducteurs du reseau d'imagerie 5 
6mettent des signaux acoustiques ei(t), ces signaux 
g6n6rent au niveau des transducteurs r du reseau cible 6 
des signaux acoustiques fr(t) s'exprimant comme suit : 

/r(/)=Jhri(t)®ei(t), 

ou ® represente l*op6rateur de convolution temporelle. 

Apr6s transfontiee de Fourier, cette Equation 
devient : 

F(fflk) = H((Dk) .E(a>k) , o£i : 

- H(cok) est la matrice de transfert, de taille 
in*n, entre les transducteurs Ti du reseau d' imagerie et les 
transducteurs Tr du r6seau cible : les composantes Hri ((Dk) 
de cette matrice sont les composantes des transform6es de 
Fourier des reponses impulsionnelles hri(t) a la frequence 
(Dk, 

- E (G>k) est un vecteur dont les composantes 
Ei(a)k) sont les composantes de la transformee de Fourier 
des signaux ei(t) susmentionnes a la frequence cok. 
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et F(a>k) est un vecteur dont les composantes 
Fj(a»k) sont les composantes de la transformfee de Fourrier 
des signaux f,(t) susmentionnes ^ la frequence cok. 

Par inversion de chaque matrice de transfert H«ak), 
5 on peut done determiner le vecteur E(o>k,j) qui est adapt! 
pour g6n6rer au niveau du point de focalisation r 
correspondant au transducteur T'r du reseau cible, un 
vecteur F(cok,j) dont toutes les composantes se rapprochent 
au mieux de I'objectif initialement fix6 (de pr6f France 
10 toutes 6gales ^ z6ro, sauf la composante d'indice j 
correspondant au transducteur T'j, qui est 6gale d 1 
lorsqu'on veut 6mettre una impulsion acoustique au niveau 
du point de focalisation r) , grace ^ la relation : 

E(ci>k, j)= H-^(fi)k) .F(Q)k, j), 
15 otl H-^ok) est la matrice inverse de H{mk) . 

H-Ucok) peut etre calculee par exemple par 
decomposition en valeurs singuli^res, ceci permettant une 
r6gularisation de 1' inversion de la matrice H (a,k) . 

Puis on determine, par transform^e de Fourier 
inverse des diff^rentes composantes Ei(fi>k,j) du vecteur 
E(cok,j), les differents signaux de reference ei(t,j) qui 
lorsqu'ils sont 6mis par les differents transducteurs Ti du 
reseau d'imagerie 5, sont adaptes pour focaliser une 
impulsion acoustique (ou le cas ech^ant un autre signal 
acoustique) au niveau du point de focalisation r. On a done 
realise une focalisation du reseau d' imagerie 5 sur chaque 
transducteur du reseau cible 6, par filtre inverse spatio- 
temporel. 

L' unite centrale CPU suit ensuite un processus 
d'apprentissage des aberrations dues ^ la paroi du crSne 3 
au niveau du reseau cible. 

Au cours de ce processus, ces aberrations sont 

consid^rees comme un f iltrt* A r•or.^r^«^ • i . 

J. -Litre a reponse impulsionnelle finie. 

Dans le domaine de Fourier, ce filtre est d6fini A 
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chaque frequence (ok une amplitude Gr((ok) et une phase 
<|)r ((Ok) . 

Pour calculer ces coefficients, on compare la phase 
et 1' amplitude de 1' ensemble des vecteurs Er. A cet effet, 
5 on commence par eliminer les dephasages introduits par les 
differences de marche entre les transducteurs d' imagerie Ti 
et les different s points de focalisation indices r. Ceci 
revient S choisir un point de focalisation particulier rO^ 
et introduire pour les autres un dephasage lineaire avec la 

10 pulsation : exp (- j Ax (rO, r, i) co) avec At (rO, r, i) = (d (rO, i) - 
d(r,i))/c ou d(r,i) est la distance entre le transducteur i 
et le point de focalisation r et c est la vitesse moyenne 
des ondes acoustiques dans le milieu a imager, en 
1' occurrence le cerveau 2. 

15 Cette correction effectuee, les differences 

d' amplitude et de phase entre les vecteurs Er sont 
attributes ^ la couche aberratrice 3 situ6es centre le 
r6seau cible. On calcule alors pour chaque point de 
focalisation r, le facteur de gain Gj (a>k) et le facteur de 

20 phase (cok) : 



Gr(fi>A) 



J 



i-i 



ou Ei* est la valeur complexe conjuguee de Ei . 

25 Ces couples {Gj (cok) , <|)j (©k) } correspondent au 

facteur d' attenuation et au dephasage relatifs introduits a 
chaque frequence par la portion de couche aberratrice 3 
situee centre le point de focalisation r. lis caracttrisent 
done finalement les aberrations introduites par la portion 

30 de couche aberratrice situee centre le reseau cible. 

On ^limine ensuite les aberrations introduites par 
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la couche aberratrice 3 situ6e centre les transducteurs 
cibles, dans 1' ensemble des p matrices H (eok) = (H ji (©k) ] , 
d6finies ci-dessus. 

Pour cela, on calcule un nouvel ensemble de 
matrices de transfert H' (cok) = [H' ji (eok) ] caractferisant la 
propagation entre le r6seau d' imagerie et le r6seau cible 
dans un milieu virtuel pour lequel seules subsistent les 
aberrations situ^es centre le r^seau d' imagerie : 

On determine alors pour chaque transducteur r du 
r6seau cible, n composantes E'i((nk,r), i variant entre 1 
et n, telles que F (cok, r) =H' (Q)k) . E' (tok, r ) , ou E'(tok,r) = 
[E'i(o)k,r)] est un vecteur h n composantes, F(o>k,r) est un 
vecteur ^ m composantes Fl(<Bk,r), 1 variant entre 1 et m, 
ces m composantes Fl (oak, r) correspondant k une 
focalisation souhait^e des ondes acoustiques k la frequence 
ok sur le point de focalisation r correspondant ^ un 
transducteur du reseau cible. 

On en d6duit, pour chaque point de focalisation r 
correspondant k un transducteur du r6seau cible, un vecteur 
de n signaux temporels de reference e' ( t, r) = (e' i (t, r) ] , i 

variant entre 1 et n, ou eV(^^)=^£'i(*wfe,r)^^'^' en notation 
complexe . 

Ces signaux de r6f6rence e'i(t,j) sont adapt6s pour 
que 1' Emission de ces derniers respectivement par les 
diff6rents transducteurs i du r6seau d' imagerie g6n6re une 
impulsion acoustique focalis^e sur le transducteur j du 
r6seau cible en 1' absence de la couche aberratrice situ6e 
centre le reseau cible. 

On notera qu'en variante, les signaux de reference 
pourraient etre determines de la fagon suivante apr6s la 
determination des reponses impulsionnelles hri(t) et des p 
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matrices de transfert H(a>k) : 

- on corrige les matrices de transfert H(a)k) pour 
s'affranchir des aberrations gen6rees par la parol 
aberratrice 3 au voisinage de chaque point de focalisation 
5 r, cette correction etant effectu6e a partir des reponses 
impulsionnelles hri(t) pr§cedemment determinees, et on 
obtient ainsi des matrices de transfert corrigees H' (a)k) , 

on determine par inversion des matrices H' (fflk) , 
pour chaque point de focalisation r correspondant a un 
10 transducteur du reseau cible, n composantes E'i((ok,r)/ i 
variant entre 1 et n, telles que F (cok, r ) -H' (G>k) . E' (cok, r ) , 
oil E' {cok,r) = [E'i(cok,r)] est un vecteur a n composantes, 
F(a)k, r) est un vecteur a m composantes FKcok^r)/ 1 variant 
entre 1 et m, ces m composantes Fl(cok,r) correspondant a 
15 une focalisation souhaitee des ondes acoustiques A la 
frequence oak sur le point de focalisation r correspondant 
a un transducteur du r6seau cible, 

et on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du r6seau 
20 cible, un vecteur de n signaux temporels e' (t,r) = 

[e'i{t,r)], i variant entre 1 et n, ou e\it,r)=^Eiia)k,r),^^'^' en 



30 



A-l 



notation complexe, les signaux e'i(t,r) etant lesdits 
signaux de reference. 

Avantageusement, lors du calcul des matrices 
25 H' (o)k) , on assimile la parol aberratrice, au voisinage de 
chaque point de focalisation r correspondant a un 
transducteur du reseau cible, a un f litre ^ reponse 
impulsionnelle finie, defini a chaque frequence cak par une 
amplitude Gr{a)k) et une phase <|>r(G)k) calculus comme suit : 



'ri{a)k).Hri{a>k) 
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H*ri d6signe la valeur complexe conjugu^e de Hri,. 

et AT(rO,r,i)=(d(rO,i)-d(r,i))/c, d(r,i) 6tant la 
distance entre le transducteur i et le point de 
focalisation r, et d(rO,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO. 

On calcule ensuite p matrices de transfert 
corrig6es H' (oak) = [H' ji (<ok) ] , oCi 

^;/(^*) = H^y(«>*).^-^«->*'^<'"'>, qui servent A determiner les 

vecteurs E'i(a)k) comme explicite precedemment et done les 

diff6rents signaux de reference e'i(t,r), r allant de 1 ^ 
m. 

L'unit6 centrale CPU precede ensuite A un 
apprentissage de focalisation en un nombre R de points de 
focalisation pr6d6termin6s situes dans le cerveau 2, 
d- indices q compris entre m+1 et m+R, cette sous -6t ape 
consistant k determiner pour chacun de ces points de 
focalisation q, en s'^loignant pas a pas des transducteurs 
du r6seau cible, des signaux de reference e'i(t,q) ^ faire 
emettre par les diff6rents transducteurs du r6seau 
d'imagerie pour gen6rer une impulsion focalis^e sur ledit 
point de focalisation q. 

Les signaux de reference e'i(t,q) sont initialement 
determines, pour chaque nouveau point de focalisation q, 
sous la forme e'i(t,q) - e' i (t+Gi (q) , qO) pour chaque point 
de focalisation q, qO etant I'indice d • un point de 
focalisation proche du- point q pour lequel le signal de 
reference a dej^ ete determine, le retard Gi (q) etant 
initialement egal A une valeur 5i (q) /c, ou c est la valeur 
moyenne de ceierite des ondes acoustiques dans le milieu, 
et 5i(q) est egal A une difference entre d' une part, une 
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distance entre le transducteur i du r6seau d' imagerie et le 
point de focalisation qO, et d' autre part, une distance 
entre le transducteur i du reseau d' imagerie et le point de 
focalisation q. 

5 Dans le cas oil certains transducteurs d' indice v du 

r6seau d' imagerie ne sont pas situes centre la couche 
aberratrice^ il est en outre souhaitable de corriger les 
signaux de r6f6rence les signaux e'v(t,q) correspondants 
par repropagation numerique des transducteurs virtuels 
10 (situes contre la couche aberratrice 3) aux transducteurs 
reels (s6par6s de ladite couche 3 par du gel 7 ou autre) ^ 
de fagon connue en soi^ par le precede inverse de celui 
decrit ci-dessus en ce qui concerne les reponses 
impulsionnelles. 

15 On fait ensuite 6mettre par les transducteurs du 

r6seau d' imagerie, par iterations, les estimations obtenues 
des signaux de reference e'i(t,q), puis on capte avec les 
memes transducteurs les signaux, Si(t,q), retrodif f us6s par 
le milieu h6t6rogene dissipatif . 

20 Puis, on modifie pour 1' iteration suivante ces 

signaux de reference e'i(t,q) de la maniere suivante : 

ou les valeurs Ti(q) et ai(q) sont un retard et un facteur 
d' amplitude correctifs, calculees pour maximiser un critere 
25 de coherence C entre lesdits signaux retrodif fuses . 

Avantageusement, le critere C pourra etre le 

suivant 



< 


n 

1=1 


2 

> 


h 


t 

=1 





, ou : 



^l(^9) = '5f(0®e/(^^) / ® representant 1' operation de 
30 convolution, 

et <> represente une moyenne temporelle. 
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Dans ce processus d' optimisation, les valeurs Ti(q) 
peuvent Stre calcul6es A chaque iteration de fagon a 
maximiser une fonction d' intercorr61ation, pour des 
transducteurs voisins du r^seau d'imagerie, des signaux 
gi(t,q) et gi+i(t,q) susmentionn^s, et les valeurs ai(q) 
peuvent Stre calcul6es de manifere ^ 6galiser sur I'indice i 
I'amplitude maximale des fonctions gi(t,q). 

On optimise ainsi les signaux de r^f^rence 
e'i(t,q), i allant de 1 ^ n, de fagon qu' ils produisent un 
signal acoustique focalis6 de fagon precise sur le point de 
focalisation q situe dans le cerveau. Ce processus 
d' optimisation a d6ja €t6 explicits plus en detail par 
Mallart et al. (The Van Cittert-Zernike theorem in pulse 
echo measurements, j. Acoust. Soc. Am. 90(5), novembre 
1991, p. 2716-2727 ; Adaptative focusing in scattering 
media through sound speed inhomogeneities : the Van Cittert 
Zernike approach and focusing criterion, J. Acoust. Soc. 
Am. 96(6), d6cembre 1994, p. 3721-3732). 

Lorsque cette optimisation est termin^e pour un 
point de focalisation q, par exemple apr^s 2 ou 3 
iterations lorsque le crit^re C a atteint une valeur 
pr6d6termin6e (notamment proche de 2/3), on passe au point 
de focalisation q+1 suivant, etc. 

Les signaux de r6f6rence e'i(t,q) ainsi obtenus 
sont m6moris6s, par exemple dans les m^moires Ml-Mn. 

Une fois l'6tape d' apprentissage termin6e, on peut 
notamment r6aliser des images 6chographiques du cerveau 2, 
6ventuellement A un rythme 61ev6 pouvant atteindre la 
cadence d'un echographe standard, par exemple 20 S 30 
images par seconde. Pour r^aliser chacune de ces images, on 
procdde comme suit, pour chaque point de focalisation q 
appartenant au champ ^ observer : 

on fait 6mettre respect ivement par les 
transducteurs Ti du r6seau d'imagerie, lesdits signaux de 
reference ei(t,q). 
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- puis on fait capter par lesdits transducteurs du 
reseau d'imagerie des signaux si(t) retrodif f us6s par le 
milieu viscoelastique, 

on convolue le signal de reference ei(t,q) emis 
par chaque transducteur du reseau d'imagerie avec le signal 
retrodiffus6 si(t) capt6 par ce transducteur, 

puis on somme les produits de convolution ainsi 

obtenus , 

On notera que les differentes operations 
susmentionn6es r6alis6es au cours de l'6tape 
d' apprentissage ou de I'etape d'imagerie peuvent Stre soit 
programm6es dans l'unit6 centrale CPU, soit realisees en 
tout ou partie par des circuits specialises. 

Par ailleurs, on notera egalement que tous les 
transducteurs Ti, T'r pourraient etre utilis6s pour 
r6aliser les images 6chographiques du cerveau. Dans ce cas, 
le r6seau d'imagerie serait confondu avec le r6seau cible 
et chacun de ces deux reseaux comprendrait tous les 
transducteurs, le f onctionnement decrit ci-dessus 
s'appliquant alors mutatis mutandis. 
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REVENDICATIONS 



1. Proced6 non invasif de focalisation d'ondes 
acoustiques dans un milieu h^tdrogene dissipatif (2,3) 
comprenant un milieu sensiblement homog^ne (2) entoure au 
moins partiellement par une couche aberra trice dissipative 
(3) qui g€nkre des aberrations dans la propagation des 
ondes acoustiques, les ondes acoustiques 6tant emises 
depuis l'ext6rieur de la couche aberratrice (3) et 
focalis6es dans le milieu sensiblement homog^ne (2) , 
earact^risd on c« qu'il comporte les Stapes suivantes : 

(a) une etape initiale de positionnement au cours 
de laquelle on fixe un nombre t sup^rieur k 2 de 
transducteurs acoustiques (Tl-Tn, T'1-T'm) dans des 
positions pr6d6termin6es k l'ext6rieur de la couche 
aberratrice (3), ces transducteurs etant en contact avec 
ladite couche aberratrice et formant au moins : 

- un r6seau d' imagerie (Tl-Tn) qui regroupe un 
nombre n compris entre 1 et t desdits transducteuts, 

- et un r6seau cible (T'1-T'm) qui regroupe un 
nombre m compris entre 1 et t desdits transducteurs, 

(b) une 6tape d' apprentissage comprenant elle-m€me 
les sous 6tapes suivantes : 

(bl) une sous-6tape d' apprentissage de 
focalisation du reseau d' imagerie sur le r6seau cible, 
sous-6tape au cours de laquelle : 

(bll) on determine des r6ponses 

impulsionnelles hri(t) du milieu heterogene dissipatif, 
respectivement entre chaque transducteur i du reseau 
d' imagerie et plusieurs points de focalisation r situ6s sur 
la couche aberratrice (3) en correspondance respective avec 
des transducteurs du reseau cible, ces reponses 
impulsionnelles etant memorisees sous forme numerique avec 
un certain echantillonnage temporel qui determine un nombre 
p de composantes f r^quentielles de la r^ponse 
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impulsionnelle, de frequences respect ives cok, i etant un 
indice compris entre 1 et n qui designe un transducteur du 
reseau d'imagerie, r etant un indice compris entre 1 et m 
qui designe un point de focalisation correspondant ^ un 
5 transducteur du reseau cible et k 6tant un indice compris 
entre 1 et p qui designs une composante f requentielle, 

(bl2) k partir de ces r6ponses 
impulsionnelles, on calcule, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 

10 cible, un ensemble de n signaux de reference temporels 
e'i(t,r), i variant entre 1 et n, tels que, si la parol 
aberratrice 6tait enlev6e au voisinage du point de 
focalisation r, 1' Emission de ces signaux de r6f6rence par 
les differents transducteurs i du reseau d' imagerie 

15 gen^rerait un signal predetermine focalis6 sur le point de 
focalisation r, 

(b2) une sous-etape de focalisation en un noxtibre 
R de points de focalisation predetermines situ6s dans le 
milieu sensiblement homogene, d' indices q compris entre m+1 

20 et m+R, cette sous-6tape consistant a determiner pour 
chacun de ces points de focalisation q, en s'eloignant pas 
^ pas des points de focalisation lam correspondant aux 
transducteurs du reseau cible, des signaux de reference 
e'i(t,q) ^ faire 6mettre par les differents transducteurs i 

25 du reseau d' imagerie pour generer une impulsion focalis6e 
sur ledit point de focalisation q, les signaux de reference 
e'i(t,q) etant determines pour chaque point de focalisation 
q en procedant comme suit : 

(b21) une premiere estimation de e'i(t,q), 

30 pour i allant de 1 ^ n, est calculee a partir d'au moins un 
signal de reference e'i(t,qO), qO etant 1 ' indice d'au moins 
un point de focalisation proche du point de focalisation q 
pour lequel le signal de reference a deja ete determine, ce 
calcul etant fait en utilisant une ceierite moyenne des 

35 ondes acoustiques dans le milieu sensiblement homogene (2), 
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(b22) on fait emettre par les transducteurs 
du reseau d'imagerie, par iterations, les estimations 
precedemment obtenues des signaux de r6f6rence e'i(t,q), 
puis on capte avec les memes transducteurs des signaux 
Si(t,q) r§trodif fuses par le milieu heterogene dissipatif, 
puis on modifie pour 1' iteration suivante ces signaux de 
reference e'i(t,q) de la mani6re suivante : 

e\it) €Ciiq)^eiit^Tiiq)) 
oh les valeurs ai(q) et Ti{q) sent un facteur d' amplitude 
et un retard correctifs, calculus pour maximiser un critere 
de coherence C entre lesdits signaux retrodif fuses, 
lesdites iterations etant arretees lorsque le critere C 
atteint un seuil predetermine, 

(b3) on memorise les signaux de reference 
e'i(t,q), au moins pour q compris entre m+l et m+R, 

(c) et une etape de focalisation au cours de 
laquelle, pour au moins un desdits points de focalisation 
q, on fait emettre respectivement par les transducteurs du 
reseau d'imagerie, lesdits signaux de reference e'i{t,q), i 
etant un indice compris entre 1 et n designant un 
transducteur du reseau d'imagerie. 

2. Procede selon la revendication 1, dans lequel au 
cours de la sous-etape (bll) , lorsqu' au moins certains 
transducteurs (Tl-Tm, T'1-T'm) sont en contact avec un 
milieu homogene intermediaire lui-meme en contact avec la 
couche aberratrice, on corrige les r6ponses impulsionnelles 
hri{t) par repropagation numerique pour simuler des 
transducteurs situes directement au contact de la couche 
aberratrice. 

3. Precede selon la revendication 1 ou la 
revendication 2, dans lequel la sous-etape (bl2) coraporte 
elle-mSme les sous-etapes suivantes : 

(bl21) on determine p matrices de transfert 
H(cok) = [Hri«i)k)], i allant de 1 a n et r allant de 1 k m. 
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ou Hri((Dk) est la valeur, a la frequence a>k, de la 
transform^e de Fourier de la reponse impulsionnelle hri(t), 

(bl22) on determine pour chaque point de 
focalisation r correspondant ^ un transducteur du r^seau 
cible, n composantes Ei(a)k, r), i variant entre 1 et n, 
telles que F (ok, r ) =H (cok) . E (ok, r ) , oii E(cok,r) = [Ei(G)k,r)] 
est un vecteur a n composantes, F(tok,r) est un vecteur ^ m 
composantes Fl(cok,r), 1 variant entre 1 et m, ces m 
composantes Fl(G)k,r) correspondant A une focalisation 
souhaitee des ondes acoustiques a la frequence cok sur le 
point de focalisation r correspondant ^ un transducteur du 
r^seau cible, 

(bl23) on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du r6seau 
cible, un vecteur de n signaux temporels e ( t , r) = [ei ( t , r ) ] , 

i variant entre 1 et n, ou ei{t,r)='^Eiicok,r).^''*' en notation 

complexe, ces signaux ei(t,r) etant adaptes pour que 
1' emission de ces derniers respectivement par les 
diff brents transducteurs i du reseau d' imagerie genere une 
impulsion acoustique focalis^e sur le point de focalisation 
r du reseau cible, 

(bl24) une sous-etape de correction des 

aberrations gen^rees par la couche aberratrice entre le 
milieu sensibleraent homogene et chaque transducteur cible 
r, ces aberrations 6tant estimees sur la base des mesures 
pr6c6demment effectuees, les aberrations ainsi estimees 
etant utilis^es pour calculer iesdits signaux temporels de 
reference e'i(t,r). 

4. Precede selon la revendication 3, dans lequel au 
cours de la sous-etape (bl22) on calcule p matrices 
H-^(o>k), respectivement par regularisation et inversion des 
matrices de transfert H{o)k), et pour chaque transducteur r 
du reseau cible, on calcule le vecteur E{a)k,r) par la 
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formule : 

E(CDk, r)=H'^(a>k) .F(a)k, j) . 

5. Proc6de selon la revendication 3 ou la 
revendication 4, dans lequel au cours de I'etape (bl22) , 
les composantes Fl(a)k, r) du vecteur F(a)k,r) correspondant 
^ la repartition spatiale du champ desir6 a la frequence 
(Ok, sont egales a 0 pour l^^r et egale ^ 1 pour l=r. 

6. Proced6 selon I'une quelconque des 
revendications 3 a 5, dans lequel au cours de la sous-6tape 
(bl24), on assimile la parol aberratrice, au voisinage de 
chaque point de focalisation r correspondant a un 
transducteur du reseau cible, a un filtre k reponse 
impulsionnelle finie, defini a chaque frequence ok par une 
amplitude Gr (ook) et une phase (t»r(CDk), la sous-6tape (bl24) 
comportant elle-meme les sous-etapes suivantes : 

(bl241) on calcule, pour chaque frequence 
cok, 1' amplitude Gr(a>k) et la phase (t)r(cok), a partir soit 
des signaux ei(t,r), soit des vecteurs E{cok, r) , 

{bl242) on calcule p matrices de transfert 
corrigees H' (tt>k) = [H' j i (cok) ] , ou 

(bl243) on determine pour chaque 
transducteur r du reseau cible, n composantes E'i(ok,r), i 
variant entre 1 et n, telles que F (cok, r) =H' (©k) . E' (cok, r) , 
oCi E' ((Dk,r)= [E'i{(Dk,r)] est un vecteur k n composantes, 
F(a)k,r) est un vecteur k m composantes Fl((ok,r), 1 variant 
entre 1 et m, ces m composantes Fl(cak,r) correspondant a 
une focalisation souhaitee des ondes acoustiques A la 
frequence cok sur le point de focalisation r correspondant 
a un transducteur du reseau cible, 

(bl244) on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du r6seau 
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cible, un vecteur de n signaux temporels de reference 
e' (t,r) = [e'i{t,r) ] , i variant entre 1 et n, oCi 

eXt^)=^Ei{cDk,r).0'^' en notation complexe. 

7- Proc6de selon la revendication 6, dans lequel au 
cours de la sous-6tape (bl241), on calcule 1' amplitude 
Gr(a)k) et la phase (|>r (tok) conune suit : 



Gr(a;*)=- 



n 

5]£/(fl>*,r0)^/*(<B»,/O) 



?i^(^0=il;(arg(£<d,*.H)))-arg(£.<a>*./-)g-^"'*"'''^)) 

10 ou : 

Ei* est la valeur complexe conjuguee de Ei, 
et AT(rO,r,i)=(d(rO,i)-d(r,i) )/c, d(r,i) etant la 
distance entre le transducteur i et le point de 
focalisation r, et d(rO,i) 6tant la distance entre le 
15 transducteur i et un point de focalisation particulier rO. 

8. Proc6d§ selon la revendication 1 ou la 
revendication 2, dans lequel la sous-etape (bl2) comporte 
elle-mSme les sous-6tapes suivantes : 

(bl21) on determine p matrices de transfert 
20 H(a)k) = [Hri(a)k) i allant de 1 a n et r allant de 1 a 

ou Hri(cok) est la valeur, a la frequence cok, de la 
transform6e de Fourier de la reponse impulsionnelle hri(t), 

(bl22' ) on corrige les matrices de transfert 
H(a)k) pour s' af f ranchir des aberrations g^nerees par la 
2 5 parol aberratrice au voisinage de chaque point de 
focalisation r, cette correction etant effectuee a partir 
des reponses impulsionnelles hri(t) precederament 
determin6es, et on obtient ainsi des matrices de transfert 
corrig^es H' (©k) , 



1 I 

2815717 



10 



30 

{bl23' ) on determine pour chaque point de 
f ocalisation r correspondant a un transducteur du r6seau 
cible, n composantes E'i(<ok,r), i variant entre 1 et n, 
telles que F(a)k,r)=H' (ok) .E' (cok.r) , ou E' ((Dk,r) = 
[E'i(G)k,r)] est un vecteur a n composantes, F{(Ok,r) est un 
vecteur a m composantes Fl(o)k,r), 1 variant entre 1 et m, 
ces m composantes Fl(cok,r) correspondant a une 
f ocalisation souhaitee des ondes acoustiques a la frequence 
(Ok sur le point de focalisation r correspondant ^ un 
transducteur du r6seau cible, 

(bl24') on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du r6seau 
cible, un vecteur de n signaux temporals 
e' (t,r) = [e'i (t,r) ] , i variant entre 1 et n, ou 

15 eV(^r)=]£]£'z(fi>^,r).g^'^' en notation complexe, les signaux 

e'i(t,r) etant lesdits signaux de r6f6rence. 

9. Precede selon la revendication 8, dans lequel au 
cours de la sous-etape (bl23') on calcule p matrices H'" 
^(©k), respectivement par regularisation et inversion des 

20 matrices de transfert H' (ok) , et pour chaque transducteur 
r du reseau cible, on calcule le vecteur E' (a)k,r) par la 
formule : 

E' (a)k,r)=H'"^(cok) .F(<Dk, j) . 

10. Proc6d§ selon la revendication 8 ou la 
25 revendication 9, dans lequel au cours de 1 ' etape (bl23'), 

les composantes Fl((Dk,r) du vecteur F(G>k,r) correspondant 
k la repartition spatiale du champ desire a la frequence 
(Ok, sont egales a 0 pour l^tr et egale a 1 pour l=r. 

11. Proc6d6 selon I'une quelconque des 
revendications 8 k 10, dans lequel au cours de la sous^ 
etape (bl22' ) , on assimile la parol aberratrice, au 
voisinage de chaque point de focalisation r correspondant a 
un transducteur du r6seau cible, k un f litre k r6ponse 
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impulsionnelle finie, defini a chaque frequence o)k par une 
amplitude Gr(cok) et une phase <|)r(o>k), la sous-etape (bl22' ) 
comportant elle-meme les sous-§tapes suivantes : 

(bl22'l) on calcule, pour chaque frequence 
ok, 1' amplitude Gr(CDk) et la phase <t)r(cok), a partir des 
r6ponses impulsionnelles precedemment determinees, 

(bl22'2) on calcule p matrices de transfert 
corrigees H' (ok) = [H' j i (cok) ] , ou 

12. Proc6de selon la revendication 11, dans lequel 

au cours de la sous-etape (bl22'l}on calcule, pour chaque 

frequence ok, 1' amplitude Gr(ok) et la phase <|>r(ok), de la 
maniere suivante : 



4^i<^^>\t.^^ ou : 

H*ri d6signe la valeur complexe conjuguee de Hri, 
et Ax(rO,r,i) = (d(rO,i)-d(r,i) )/c, d(r,i) etant la 
distance entre le transducteur i et le point de 
focalisation r, et d(rO,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO. 

13. Procede selon 1 ' une quelconque des 
revendications prec6dentes, dans lequel au cours de I'etape 
(c) , on fait suivre la sous-etape (cl) par les sous-etapes 
suivantes : 

(c2) on fait capter par lesdits t ransducteurs du 
reseau d'imagerie des signaux Si(t) retrodif fuses par le 
milieu heterogene dissipatif. 
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(c3) on convolue le signal de reference 6mis par 
chaque transducteur du reseau d'imagerie avec le signal 
retrodiffus§ capt6 par ce transducteur, 

(c4) puis on somme les produits de convolution 
5 ainsi obtenus, 

l'6tape (c) 6tant renouvel6e pour une plurality de points 
situ6s dans le milieu sensiblement homogfene. 

14. Proc6d4 selon I'une quelconque des 
revendications pr6c6dentes, dans lequel au cours de la 

10 sous-6tape (b21) , la premiere estimation de chaque signal 
de r6f6rence est e^i(t,q) = e' i ( ts+Gi (q) , qO ) pour chaque 
point de focalisation q, qO etant I'indice d'un point de 
focalisation proche du point q pour lequel le signal de 
reference a d6jA et§ determine, 9i(q) etant un retard egal 

15 a une valeur 5i (q) /c, ou c est la c61erit6 moyenne des 
ondes acoustiques dans le milieu, et 5i(q) est 6gal k une 
difference entre d'une part, une distance entre le 
transducteur i du reseau d'imagerie et le point de 
focalisation qO, et d' autre part, une distance entre le 

20 transducteur i du r6seau d' imagerie et le point de 
focalisation q, 

15. Proced6 selon I'une quelconque des 
revendications pr6c6dentes, dans lequel au cours de la 
sous-etape (b2), lorsqu'au moins certains transducteurs 

25 d'indice v du reseau d' imagerie ne sont pas directement au 
contact de la couche aberratrice, on corrige les signaux 
e'v{t,q) * correspondants par repropagation numerique pour 
simuler des transducteurs places en contact direct avec la 
couche aberratrice. 

16. Precede selon I'une quelconque des 
• revendications pr6c6dentes, dans lequel au cours de la 

sous-6tape (b22), on recherche les valeurs ai(q) et Xi (q) 
pour maximiser le critere de coherence C suivant : 
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n 



ou 



i=i 

g^,(/,9) = 5;(0®^,(/,g) / ® representant 1' operation de 
convolution, 

et <> represente une moyenne temporelle. 
5 17. Proc6d6 selon la revendication 16, dans lequel 

au cours de la sous-etape {h22) , les valeurs Ti(q) sont 
calcul6es en maximisant une fonction d' intercorr61ation, 
pour des transducteurs voisins du r6seau d' imageries des 
signaux gi(t,q) et gi+i(t,q). 
10 18. Proced6 selon la revendication 16 ou la 

revendication 11, dans lequel au cours de la sous-6tape 
(h22) , les valeurs ai(q) sont calculees de maniSre a 
egaliser sur I'indice i 1' amplitude maximale des fonctions 
gi(t,q) . 

15 19. Proc6de selon l*une quelconque des 

revendications pr6cedentes, dans lequel la sous-etape (b22) 
relative S chaque point de focalisation q est realis6e 
immfediatement apres la sous-etape (b21) relative au m6me 
point de focalisation q. 

20 20. Proced6 selon I'une quelconque des 

revendications pr6cedentes, dans lequel le milieu 
heterogene dissipatif est constitue par le cerveau entoure 
par le crSne. 

21. Proced6 selon I'une quelconque des 
25 revendications pr6c6dentes, dans lequel : 

- soit le reseau d'imagerie et le reseau cible 
sont deux reseaux distincts disposes de part et d' autre du 
milieu heterogene dissipatif^ 

- soit tous les transducteurs appartiennent a la 
30 fois au reseau d'imagerie et au reseau cible. 

22. Proc6de selon I'une quelconque des 
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revendications precedentes, dans lequel les ondes 
acoustiques sont des ondes ultrasonores . 

23. Dispositif (1) adapts pour mettre en oeuvre un 
proc6d6 selon I'une quelconque des revendications 
pr6c§dentes, ce dispositif comportant un norabre t sup6rieur 
a 2 de transducteurs acoustiques (Tl-Tn, T'1-T'm) destines 
^ Stre fix6s dans des positions pr^determin^es k 
l'ext6rieur de ia couche aberratrice (3), ces transducteurs 
6tant commandos par au raoins une unit6 centrale 
61ectronique (CPU) et formant au moins : 

- un r6seau d' imagerie (Tl-Tn) qui regroupe un 
nombre n compris entre 1 et t desdits transducteurs, 

- et un r^seau cible (T'1-T'm) qui regroupe un 
nombre m compris entre 1 et t desdits transducteurs, 
l'unit6 centrale 61ectronique 6tant adapt6e pour suivre les 
Stapes suivantes : 

(b)une 6tape d' apprentissage comprenant elle-mSme 
les sous Stapes suivantes : 

(bl) une sous-etape d' apprentissage de 
20 focalisation du reseau d' imagerie sur le r6seau cible, 
6tape au cours de laquelle : 

(bll) on determine des r6ponses 

impulsionnelles hri(t) du milieu h6t6rog6ne dissipatif, 
respectivement entre chaque transducteur i du reseau 
25 d' imagerie et plusieurs points de focalisation r situ6s sur 
la couche aberratrice en correspondance respective avec des 
transducteurs du reseau cible, ces r^ponses impulsionnelles 
etant memorisees sous forme numerique avec un certain 
echantillonnage temporel qui determine un nombre p de 
composantes f r^quentielles de la r6ponse impulsionnelle, de 
frequences respectives cok, i 6tant un indice compris entre 
1 et n qui designe un transducteur du reseau d' imagerie, r 
6tant un indice compris entre 1 et m qui d6signe un point 
de focalisation correspondant k un transducteur du reseau 
cible et k 6tant un indice compris entre 1 et p qui d6signe 
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une composante f r6quentielle, 

(bl2) a partir de ces reponses 

impulsionnelles^ on calcule^ pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
5 cible, un ensemble de n signaux de reference temporels 
e'i(t,r)^ i variant entre 1 et n, tels que, si la parol 
aberratrice 6tait enlevee au voisinage du point de 
focalisation r, 1' emission de ces signaux de reference par 
les diff6rents transducteurs i du reseau d' imagerie 

10 g§n6rerait une impulsion acoustique focalisee sur le point 
de focalisation r, 

(b2) une sous-6tape de focalisation en un nombre 
R de points de focalisation predetermines situes dans le 
milieu sensiblement homog6ne, d' indices q compris entre m+1 

15 et m+R, cette sous-etape consistant k determiner pour 
chacun de ces points de focalisation q, en s'61oignant pas 
a pas des points de focalisation 1dm correspondant aux 
transducteurs du reseau cible, des signaux de reference 
e'i(trq) d faire emettre par les differents transducteurs i 

20 du r6seau d' imagerie pour generer un signal pred6teirmin6 
focalise sur ledit point de focalisation q, les signaux de 
reference e'i(t,q) 6tant determines pour chaque point de 
focalisation q en procedant comme suit : 

(b21) une premiere estimation de e'i(t/q), 

25 pour i allant de 1 k 4, est calculee a partir d'au moins un 
signal de r6f6rence e'i(t,qO), qO etant 1 ' indice d'au moins 
un point de focalisation proche du point de focalisation q 
pour lequel le signal de reference a deja 6te determine, ce 
calcul etant fait en utilisant une celerite moyenne des 

30 ondes acoustiques dans le milieu sensiblement homogene (2) , 

(b22) on fait emettre par les transducteurs 
du reseau d' imagerie, par iterations, les estimations 
precedemment obtenues des signaux de reference e'i(t,q), 
puis on capte avec les memes transducteurs des signaux, 

35 Si(t,q) retrodif fuses par le milieu heterogene dissipatif. 
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puis on modifie pour 1' iteration suivante ces signaux de 
reference e'i(t,q) de la mani6re suivante : 

oh les valeurs ai(q) et Ti(q) sent un facteur amplitude 
et un retard correctifs, calcules pour maximiser un critere 
de coherence C entre lesdits signaux retrodif f us6s, 
lesdites iterations etant arretees lorsque le critere C 
atteint un seuil predetermine^ 

(b3) on memorise les signaux de reference 
e'i(t,q), au moins pour q compris entre m+1 et m+R, 

(c) et une etape de focalisation au cours de 
laquelle, pour au moins un desdits points de focalisation 
q/ on fait 6mettre respectivement par les transducteurs du 
r6seau d'imagerie, lesdits signaux de reference e'i(t,q), i 
etant un indice compris entre 1 et n designant un 
transducteur du reseau d'imagerie. 
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